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Abstract 

In the event that R-Parity conserving supersymmetry (SUSY) is discovered at the 
LHC, a key issue which will need to be addressed will be the consistency of that signal with 
astrophysical and non-accelerator constraints on SUSY Dark Matter. This issue is studied 
for the SPA benchmark model based on measurements of end-points and thresholds in the 
invariant mass spectra of various combinations of leptons and jets. These measurements 
are used to constrain the soft SUSY breaking parameters at the electroweak scale in a 
general MSSM model. Based on these constraints, we assess the accuracy with which the 
Dark Matter relic density can be measured. 
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1 Introduction 


The complementarity of LHC SUSY measurements and direct and indirect searches for neu- 
tralino Dark Matter is an important topic to study given the increasingly strong astrophysical 
evidence for cold Dark Matter in the universe PHEl. Assuming that R-Parity conserving 
SUSY is discovered at the LHC, an interesting question will arise regarding the compatibility 
of that signal with existing relic density constraints (e.g. 0.094 < < 0.129 at 2a from 

WMAP data and the implications it has for terrestrial Dark Matter searches. 

In a previous paper we addressed these issues within the context of the minimal su¬ 
pergravity (mSUGRA) class of SUSY models incorporating gravity-mediated SUSY breaking 
1^. There is a great advantage in studying models such as the mSUGRA ones, where a 
definite pattern of SUSY breaking is imposed at the unification scale. These models are in 
fact described with only a small number of independent parameters. R is therefore possible, 
through a limited number of measurements to fully constrain the model, and extract very 
precise predictions for physical quantities related to Dark Matter. However, the mSUGRA 
model is strongly constrained by the WMAP data, and we do not know whether from the LHC 
measurements an unambiguous determination of the high scale behaviour of the SUSY theory 
will be possible. It is therefore interesting to study a generic model in which the weak-scale 
SUSY breaking parameters are independent, such as the Minimal Supersymmetric Standard 
Model (MSSM), and verify to what level the LHC measurements can constrain the neutralino 
relic density and the cross-sections for direct Dark Matter detection. 

Some recent papers 013, in the framework of International Linear Collider (ILC) studies, 
based on generic and sometimes restrictive assumptions on the measurement potential of 
the LHC, conclude for a very limited potential of the LHC experiments in constraining the 
neutralino relic density in a generic MSSM. The aim of this paper is to investigate in detail, by 
carefully taking into account published studies, what are the effective limitations of the LHC 
measurements in this held for a specihc representative model for which sufficient information 
is available. 

We choose to study one particular point of the MSSM parameter space, which was adopted 
as a benchmark point by the SPA group [H]- This model is dehned in terms of the parameters 
of the mSUGRA model (mo = 70 GeV, mi /2 = 250 GeV, Aq = —300 GeV, tan/3 = 10, /r > 0). 
This is a modihcation of the point SPSla, essentially achieved by lowering rriQ from 100 to 
70 GeV, originally dehned in Ref. 0, to take into account the WMAP results. It is also very 
near to Point B' of m- This model lies in the ‘bulk’ region of the mo — rni /2 mSUGRA 
plane where the relic density is reduced to a small value by annihilation to leptons or 
neutrinos via t-channel slepton, stau and sneutrino exchange. The values of the MSSM soft 
parameters and of the sparticle masses were computed with the program ISASUSY 7.71 ini, 
taking the tree-level values for the sparticle masses, and are given in Table 0 The Lightest 
Supersymmetric Particle (LSP) relic density was computed based on these inputs using the 
program micrOMEGAs 1.3.6 [13, and the resulting value is H^/i^=0.108, well within the 
WMAP range. The main annihilation processes contributing to the calculation of the relic 
density are given in Table (3 As discussed above we take this point as a benchmark, but we 
analyze it without assuming a predefined relationship among the SUSY breaking parameters. 

In Section 2 we describe the LHG measurements which can be used to constrain SUSY 
at the considered model point. In Section 3 we develop a strategy to constrain the MSSM 
parameters relevant to the calculation of the neutralino relic density. In Section 4 we discuss 
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the use of these constraints to calculate the x? relic density. Finally in Section 5 we review 
the results and discuss the general applicability of the technique. 


Sparticle 

mass (GeV) 

Sparticle 

mass (GeV) 

Xi 

97.2 

X2 

180.1 

X3 

398.4 

X4 

413.8 

h 

189.4 


124.1 

Tl 

107.7 

T2 

194.2 

h 

347.3 

h 

562.3 

UL 

533.3 

9 

607.0 

h 

116.8 

A 

424.6 


Table 1: Masses of the sparticles in the considered model as calculated at tree level with 
IS A JET 7.71 m 


Process 

Fraction 

XiXi ^ 

40% 

XiXi ^ t + t - 

28% 

XiXi ^ 

3% 

X?fi ^ Zr 

4% 

X?ri ^ At 

18% 

flfi TT 

2% 


Table 2: Fractional contribution of different annihilation processes to the LSP relic density in 
the considered model. The relie density was ealculated with the program micrOMEGAs 1.3.6 

m- 


2 The measurement of SUSY parameters at the LHC 

The development of techniques for measuring parameters characterizing SUSY models has 
been the subject of much investigation in the last few years, as documented in Ref. umiiin], 
and is still a very active field of investigation. 

The basic issue is that the presence of two invisible particles in the final state renders the 
direct measurement of sparticle masses through the detection of invariant mass peaks impossi¬ 
ble. Alternative techniques have therefore been developed, based on the exclusive identification 
of long cascades of two body-decays. It was demonstrated [El El that if a sufficiently long 
chain can be identified (at least three successive two-body decays), the thresholds and end¬ 
points of the various possible invariant mass combinations among the identified products can 
be used to achieve a model-independent measurement of the masses of the involved sparticles. 
Once the masses of the lighter sparticles are obtained with this procedure, in particular the 
mass of the LSP, additional sparticle masses can be measured through the identification of 
other shorter decay chains. This program has been carried out recently for the SPSla model 
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point m, assuming the performance of the ATLAS detector. This study results in a number 
of measurements of observables which are related to the masses of the sparticles by known 
algebraic relations. For the work presented here, we assume for the corresponding observables 
in the SPA model the same errors as the ones obtained in the full analysis for Point SPSla. 
The position of the kinematic edges is quite similar in the two points, due to the fact that 
only the slepton spectrum is modified by the 30% change in mo, and the total production 
cross-section is somewhat higher in the new point, which should result is somewhat smaller 
statistical errors. We further checked that the mass spectrum for the SPA model does not 
present degeneracies in the sparticle masses which would severely reduce the transverse mo¬ 
menta of the visible decay products and thence the experimental selection efficiency. The 
considered measurements are give in Table |31 The meaning of the different observables, and 
their expression in terms of sparticle masses is given in Ref. m- The scale error is derived 
from the assumption of a control on the lepton energy scale at the level of 0 . 1 %, and of the 
jet energy scale at the level of 1 % M- Following the 1 % jet scale uncertainty reflects in 
a ~ 0.5% uncertainty on the position of the edges involving jets and leptons. 


Variable 

Value (GeV) 

E 

Stat-|-Sys (GeV) 

rrors 

Scale (GeV) 

Total 

r^max 

81.2 

0.03 

0.08 

0.09 

^max 

425.3 

1.4 

2.1 

2.5 

fyy^lOW 

^iq 

266.9 

0.9 

1.3 

1.6 

hiqh 

365.9 

1.0 

1.8 

2.1 

^min 

^Uq 

207.0 

1.6 

1.0 

1.9 

m{iL) - m{xi) 

92.3 

1.6 

0.1 

1.6 


315.8 

2.3 

0.3 

2.3 

^max 

iiLrr 

62.2 

5.0 

0.3 

5.0 


Table 3: Summary table of the SUSY measurements which can be performed at the LHC 
with the ATLAS detector m used in this analysis. The central values are calculated with 
ISASUSY 7.71, using the tree-level values for the sparticle masses. The statistical errors are 
given for the integrated luminosity of 300 fb~^. The uncertainty in the energy scale is taken 
to result in an error of 0.5% for measurements ineluding jets, and of 0.1% for purely leptonic 
mesurements. 


Some additional measurements are considered in this analysis: 

• The mass of the light Higgs boson from the decay /i —> 77 which for 100 fb“^, and for 
m/i = 114 GeV, has a statistical uncertainty of ~ 0.5 GeV m- 

• The ratio of branching ratios: 

BRixl ^ 4^)/HR(xO ^ fir). 

No detailed experimental study is available on this measurement. From a particle-level 
evaluation of the number of events involved, the error on this quantity will be domi¬ 
nated by the systematic uncertainty in the experimental efficiency for the fi channel, in 
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particular near threshold where the contribution to the visible ditan mass distribution 
is significant. We assume here a conservative error of 10%. 

• The constraints from searches for heavy Higgs bosons of MSSM, which will be discussed 
in detail in the following. 

Two different types of uncertainties are quoted in Table |21 the combined statistical and 
systematic uncertainties estimated for each measurement, and general errors on the scales of 
lepton and hadron energy measurement, which affect all the measured quantities in the same 
way. Since in many cases the scale uncertainties are dominant it is necessary to take into 
account the correlations between the different measurements when extracting the constraints 
on the model parameters. In order to take into account the correlations we use a Monte Carlo 
technique relying on the generation of simulated experiments sampling the probability density 
fnnctions of the measured observables. In frequentist statistics, confidence bands describe the 
probability that an experiment in a set of identical experiments yields a given value for the 
measured quantities. 

We proceed in the following way: 

1. An ‘experiment’ is defined as a set of measurements, each of which is generated by 
picking a value from a gaussian distribution with mean given by the central value given 
in Table 01 The correlation due to energy scale is taken into account by using a second 
gaussian distribution for the energy scale, and using the same random number for all 
the measurements sharing the same scale. 

2. For each experiment, we extract the constraints on the MSSM model as we will describe 
in the following. 

We obtain as a result of this calculation a set of MSSM models, each of which is the “best” 
estimate for a given Monte Carlo experiment of the model generating the observed measure¬ 
ment pattern. For each of these models the properties of the LSP Dark Matter candidate 
and other SUSY particles involved in Dark Matter annihilation can then be calculated, with 
the spread of obtained results being interpreted as the level of precision with which these 
properties can be measured by the LHC. As we are working in a general MSSM, some of 
the parameters will only be loosely constrained, if at all, by the measurements. The spread 
obtained from the variation of the unconstrained parameters needs to be studied in detail to 
assess the final prediction precision on DM characteristic. 

3 Extraction of MSSM parameters 

In order to extract the parameters of the MSSM we proceed in a stepwise fashion. We first 
extract the measurement of the sparticle masses from the measured edges. We do not address 
here the issue of possible ambiguities which can arise on the identity of the particles involved 
in the decay chains as discussed in m and EOl. The procedure yields an error of ~9 GeV for 
the masses of the sparticles relevant to this study. We show in the left side of Figure ^ the 
distribution of the measured x? masses for a set of Monte Carlo experiments. The error on 
the mass is basically independent on the mass of the sparticle because the errors on masses 
are strongly correlated. This is because of the algebraical structure of the formulas in m, 
which causes the absolute scale to be much more loosely constrained than the difference 
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Figure 1: Distribution of the measured value of m(x?) (left) and of the difference m(x 2 ) “ 
m(xi) (right) for a set of Monte Carlo experiments, each corresponding to an integrated statis¬ 
tics of 300 fh~^. 


between masses. In particular, compared to the absolute error of ~ 9 GeV quoted above, 
the error on the difference — m-(xi) is ~ 200 MeV , as shown on the right side of 

Figure ^ thanks to the excellent precision of the measurement of the edge position. 

The calculated precision on m(fi) — m(xi) is ~ 2.5 GeV, scaling approximately linearly with 
the assumed error on . 

After extracting the mass measurements, the next step consists of constraining the neu- 
tralino sector, and calculating the components of the neutralino mixing matrix, which is a 
necessary ingredient in all calculations on the neutralino annihilation rate. Once this is fixed, 
we go on to constrain the slepton sector. The next step is understanding the constraints in 
the Higgs sector in order to evaluate possible contributions to the annihilation rate of Higgs 
exchange diagrams. By putting all of the information together we are finally able to give 
an estimate on the predictive power of the LHG data for a specific scenario in which the 
neutralino annihilation is dominated by processes involving sleptons. 

3.1 Constraints on the neutralino mass matrix 

Based on the expected LHG measurements, the masses of three neutralinos are measured: Xi? 
X 2 , X 4 - III the MSSM, assuming a CP—conserving scenario, the mass eigenstates x? (^=lj2,3,4) 
result from the mixing of the SUSY partners of the neutral gauge and Higgs bosons: 
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through a mixing matrix defined as: 


M = 


Ml 0 

0 -M 2 

—mz cos j3sw fnz cos f5cw 
y m^sin/3svF —m^sin/Jcvi/ 


—mz cos j3sw mz sin Psw \ 
mz cos Pcw —mz sin Pcw 
0 —/r 




0 


( 1 ) 




The above matrix is built out of four MSSM parameters: Mi and M 2 are respectively the 
U(l) and SU(2) gaugino masses, /x is the Higgsino mass parameter, tan/? = V 2 /V 1 is the ratio 
of the vacuum expectation values of the two Higgs doublets of the model. The additional 
parameters, sw and cw, respectively the sine and cosine of the electroweak mixing angle 
and mz, the mass of the Z boson are precisely known from measurements in the Standard 
Model. Given that we only have three input parameters in the three masses, we miss one 
parameter to fully solve the neutralino matrix. We take this parameter to be tan/3, as the 
other parameters relate directly to the neutralino masses. We show in Figure|21the distribution 
of the reconstructed values of Zn and Z 13 , for tan/3 = 10, where the composition of the x? is 
written as: 

X? = ZnB + Zi21H3 + Z 13 B? + ^14^2° 

The spread from the experimental error is 0.02% for the bino component, and 1.5% for the 



Figure 2: Distributions of the measured values of Zn, the bino component 0 /Xi (Wl) o.nd 
of Z\ 3 , the Higgsinoi component o/x? (right) assuming a fixed value tan/3 = 10 for a set of 
Monte Carlo experiments, each corresponding to an integrated statistics of 300 fb~^. 


other components. The dependence of the four components of the LSP on the assumed value 
of tan/3 is shown in Figure 01 and for the assumed range 3-30 is 0.8% for Z\\ and ~15% for 
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0.99 

0.988 

0.986 





tanp 


Figure 3; Average values of the four eomponents of Xi as a function of the assumed value of 
tan/3. The averages are performed on sets of Monte Carlo experiments, each corresponding to 
an integrated statistics of 300 fb~^. 
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Zi 3 . A much larger variation with tan/3 is observed for the smaller components Z 12 and Z 14 . 
The sensitivity is higher for low values of tan (5. A first consideration from this measurement 
is that already at this stage the data tell us that the neutralino is an almost pure bino, with 
only a small Higgsino component. Annihilation channels involving the exchange of sfermions 
in the t-channel will dominate, unless the masses of the heavy Higgses are such that they 
resonantly enhance the relevant annihilation processes. 

In the calculation above Mi and M 2 are assumed as independent, since we are working on a 
general MSSM model. It is however customary in MSSM studies to assume that the ratio of 
the gaugino masses is equal to the ratio of the respective coupling constants, giving the relation 
M 2 = 5/3 tan^ Under this assumption we can invert the gaugino mixing matrix using 

as constraints the mass differences m{^) — m(y 5 ) m(x 4 ) — instead of the measured 

masses, and extract the values of M 2 and /x, always for fixed tan /3. Prom these one can in 
turn calculate the values of the gaugino masses. In this case the absolute scale is fixed in a 
much more precise way, and the resolution on m{xi) is ~ 250 MeV. We will keep as a baseline 
for our study the completely unconstrained model, but we will also show the results when the 
gaugino mass unihcation is assumed. 

3.2 Constraints on the slepton sector 

Once the neutralino mass is extracted, the slepton sector can be considered. No mixing is 
assumed here in the selectron and smuon case. It was shown in m that the difference in mass 
of the electron and muon should result in different mixings for the selectron and the smuon, 
resulting for the considered model point in a small mass difference and in a ~ 8 % difference 
between BR{x 2 and BR{x 2 We assume for this analysis that no mixing 

takes place in the selectron and smuon sector. 

Due to the hypercharge difference, the LSP pair annihilation cross section through t- 
channel exchange of in is 16 times bigger than that for in exchange when m{iR)n = 'nn{ii,) and 
the LSP is dominantly bino. Therefore if the chirality of the slepton observed in the X 2 decay 
chain can not be determined, this could cause a very large uncertainty in the relic density 
calculation. The issue is however discussed in EH, where it is shown that the asymmetry 
distribution proposed in m is sensitive to the chirality structure of the slepton. Those 
studies were performed for the SPSla point, so their results can be safely extended to the 
model addressed in the present study, and we can assume for the present study that the 
chirality of the lighter slepton can be determined. 

In the stau sector, the ratio BR{x 2 ^R^)/BR{x 2 Bt) is a function of the neutralino 
mixing matrix, of m(x 2 )j tan/3 and 9^, the mixing angle between f/j and fi. By 

assuming a given value for tan (3 we can therefore from this measurement extract the value of 
Or- The distribution of the measured 6^, as well as the dependence on the input tan/3 value 
are are shown in Figure01 The experimental uncertainty is ~2%, whereas 9^- varies by ~ 35% 
over the considered tan /3 range. 

The only missing parameter in order to be able to calculate the neutralino annihilation 
processes involving sleptons is the mass of the f 2 slepton. Natural bounds can be imposed to 
this mass. First of all we must have m(f 2 ) > mix^ — otherwise the r 2 contribution 

would likely be observed in the rr edge from the x^ decay. Second, from m(fi), 9r, tan/3 
and m{f 2 ) the value for the trilinear coupling A^- can be calculated. Large values of the Ar 
parameter could induce charge breaking minima (CCB) due to the vacuum expectation values 
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Figure 4: Left: distribution of the measured value of Or, the f mixing angle assuming a fixed 
value tan /3 = 10 for a set of Monte Carlo experiments, each corresponding to an integrated 
statistics of 300 fb~^. Right: average values of Or, as a function of the assumed value o/tan/3. 
The averages are performed on sets of Monte Carlo experiments, each corresponding to an 
integrated statistics of 300 fb~^. 


of the charged r scalars. Typical constraints on Ar from these considerations m would give 
an upper limit on Ar of ~500 GeV, very near to the actual Ar value for the considered model. 
The conditions on CCB minima resulting in constraints on the MSSM are in general neither 
necessary nor sufficient to give an acceptable vacuum structure to the theory, and it has been 
suggested m that the constraint from m be relaxed through an empirical multiplicative 
factor. Following these considerations we conservatively fix here a limit of 5 TeV to Ar, 
which, in the considered model, for tan/3 = 10 corresponds to a limit m{f 2 ) < 250 GeV. 

3.3 Constraints from the Higgs sector 

As discussed in the previous sections, the measurements considered up to now do not allow 
us to constrain the tan /3 parameter. The main constraints come from the Higgs sector. In 
particular, the measurement of the mass and of the production rate of one of the heavy Higgs 
would define a confidence region in the canonical m{A) — tan (3 plane. The region where heavy 
Higgses can be discovered through their decays to Standard Model particles are shown in 
Figure ini from M- 

The plot was obtained on specific assumptions on the SUSY spectrum, i.e. a SUSY mass 
scale Msusy of 1 TeV, and maximal mixing, i.e. p, <C Msusy and the trilinear couplings 
A = y/^MsuSY- These assumptions mostly affect the considerations related to the light 
Higgs, and only lightly affect the discovery region for H/A and . Given that the mixing 
in the chosen model is less than maximal {At = —450 GeV), the limit on tan/3, tan/3 > 3 
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Figure 5: Reach of the ATLAS experiment in the m(A) — tan/3 plane for an integrated lumi¬ 
nosity of 300 fb~^. For each region in the plane, the detectable Higgs bosons are marked. 


shown in Fignre|Slcan be used as lower limit of the tan/3 range considered in this analysis. 
The considered model, corresponding to = 424 GeV and tan/3 = 10 lies outside of the 
region where a heavy Higgs can be discovered, we will limit onr analysis to the complement 
of this region in the m{A) — tan /3 plane. 

A stronger constraint can be obtained by the fact that the light Higgs h will be discovered at 
the LHC and its mass measured. From this measurement a confidence band in the m{A)—tan (3 
plane can be defined. 

Finally, since the SUSY spectrum is largely known, we can investigate whether the heavy 
Higgses can be detected either a) in the cascade decay of a sparticle, or b) throngh their decay 
into a sparticle, or a bonnd can be extracted from the non-observation. 

For the cascade decays, no heavy Higgs appears in the decay chains, as all the decays of 
neutralinos/charginos into heavy Higgses are kinematically closed. We can however investigate 
up to which Higgs masses these decays would be open. The best candidates wonld be the 
decays: ^4(3) ^ ^ xf-^/H with subsequent decays of the A/H into bb. These 

decays have been studied in [ISj where it is shown that in favourable conditions a peak in bb 
distribution can be observed. The kinematic limits in the considered model are: 
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• m{AlH) < 315 GeV for ^4(3) ^ Xi^/H 

• m{A/H) < 230 GeV for X4(3) ^ xl^/H and xf ^ xf^/H. 

We evaluated the number of events in which a H/A ^ bb decay is produced in the cascade 
decay of squarks and gluinos by generating with HERWIG m a sample of events for the chosen 
model, and computing the fraction of events containing either of Xs^X^^xt^- We thereafter 
computed with IS ASUS Y m the BR of the charginos/neutralinos for the considered model 
and different values of m{A), and convoluted the results with the total SUSY cross-section 
calculated at NLO with PROSPINO [13 • The expected numbers of events as a function 
of m{A) for an integrated luminosity of 300 fb“^ are shown in Figure El The curve clearly 
shows the drop in number of events when the decays to X2 xt become kinematically 
inaccessible, and the fast decrease to zero after 300 GeV. For m{A) = 300 GeV, the number 



m(A) (GeV) 


Figure 6: Number of expected events containing a H/A ^ bb decay in the cascade of squarks 
or gluinos, as a function ofm{A). The assumed integrated luminosity is 300 fb~^ 


of events contributing to the bb peak is ~1500. It should therefore be possible to put a lower 
bound on the heavy Higgs mass of approximately 300 GeV. In order to verify whether this 
is actually possible, a detailed experimental analysis is needed, outside of the scope of the 
present study. 

For the SUSY decays to heavy Higgses, two channels have been identified as particularly 
promising, and have been the subject of detailed experimental analyses 

• ^ Xi xf ^ 
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The first channel is hopeless for the considered model, as the value of the branching ratio for 
the charged Higgs decay to three leptons through a chargino-neutralino pair is of order a few 
10 -®. 

For the channel A/H ^ X 2 X 2 ’total BR to 4 leptons through the second neutralino is 
1.36x10“'^ (0.46x10“^) respectively for the A {H). We considered the A and H production 
both through gluon fusion, and in association with two b quarks, based on the available 
NLO calculations ED EH- A total of 40 events per experiment is produced for an integrated 
luminosity of 300 fb“^. The typical analysis efficiency for this signal is 30%, with the reduction 
of the SM background to negligible levels. The dominant SUSY background can be further 
reduced with respect to previous analyses by the possibility of performing a full reconstruction 
of the A/H mass peak, as discussed in ESj- It is therefore thinkable that the heavy Higgs can 
be discovered in this channel. However, given the small number of events considered, a very 
careful background study, taking into account the details of the detector performance would 
be required. 

In conclusion, three scenarios can be envisaged: 

• A scenario in which the only constraints on the m{A) — tan/3 plane are provided by the 
measurement of the light Higgs h, and by the non-discovery of the heavy Higgses in their 
Standard Model decay modes 

• A scenario in which a lower limit of approximately 300 GeV on the H/A mass can be 
set by the non-observation of the heavy Higgses in the SUSY cascade decays. 

• A scenario in which the H/A is discovered in its X 2 X 2 decay mode. 

In the following we will evaluate the constraints on taking into account all the three 

scenarios. 

4 Calculation of the relic density 

We can at this point calculate using the micrOMEGAs program m the LSP relic density 
for each of the Monte Garlo experiments defined in Sectional The parameters derived 
from the measurements given in Table 01 used to create the Monte Garlo experiments are the 
masses of x?) X 2 ) Xi^ ^L, and fi and the stau mixing angle 9r. From these, for each 
experiment the soft MSSM parameters to give as input to micrOMEGAs are calculated. The 
nominal values are assumed for all the other MSSM parameters, in particular m{A), tan/3 
and m(f 2 ). 

The resulting distributions of the calculated are given in Figure [3 for two values of 

the assumed uncertainty on the position of the rr edge, respectively 5 GeV and 0.5 GeV. The 
error is respectively ~20% (10%) for an uncertainty of 5(0.5) GeV. This uncertainty is quoted 
in Table El as 5 GeV, and translates to an uncertainty on the mass difference m{fi) — m(x5) 
of 2.5 GeV. The measurement precision on this difference has been found in [34j to be the 
dominant factor in the determination of the relic density. No detailed study is available on 
the precision with which the LHG experiments can define the rr edge, and 5 GeV was quoted 
in m as a conservative upper limit on this figure. It is therefore interesting to treat this 
uncertainty as a parameter of the analysis. We have therefore re-evaluated the neutralino 
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relic density as a function of the error on the rr edge position. The results are shown as full 
circles in Figure |S1 The overall uncertainty depends on the rr measurement as long as the 
uncertainty is above ~ 1 GeV, which therefore should be taken as the goal for the systematic 
control to be achieved on the measurement of this variable. 

The remaining uncertainty of 10% is determined by the ~10% error on the Xi mass, as 
can be seen in Figure IHl produced for a{m{TT)) = 1 GeV, where the predicted relic density 
is shown as a function of the measured mass. If we assume gaugino mass universality, as 
discussed in Section |31 the uncertainty on the Xi mass is reduced to ~ 250 MeV. We show as 
full squares in Figure |H1 the evolution of the uncertainty as a function of the error on the rr 
edge position under this assumption. For an uncertainty on m(rr) of 1 GeV, the measurement 
error is reduced to ~ 4%. 

We have also studied the effect of changing the systematic uncertainty on the measurement 
of BR{x 2 ^r^)/BR{x 2 ^it) from 10% to 1%. The effect has been shown to be negligible 
as compared to the other uncertainty sources. 






Figure 7: Distributions of the predicted relie density incorporating the experimental 

errors. The distributions are shown for an assumed error on the rr edge respeetively of 5 GeV 
(left) and 0.5 GeV (right). 


We further studied in detail the dependence of the relic density prediction on the parame¬ 
ters which are only very loosely constrained (if at all) by the LHG data. The squark and gluino 
masses, will be measured with a precision of a few percent at the LHC, but for the masses 
considered, the relic density is essentially independent from their value. The only loophole 
could be a light stop. We have studied the dependence of the predicted relic density on the 
mass of the stop. We found that a light stop would only contribute to the relic density if its 
mass were below ~140 GeV. A detailed study on a similar model point, given in jSH], shows 
that the lighter stop, which has a mass of ~ 400 GeV, can be observed in the cascade decays 
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Figure 8: Fractional uncertainty on the predicted relic density from the experimental 

measurement, as a function of the assumed uncertainty on the position of the rr edge respec¬ 
tively for the completely unconstrained MSSM (black circles), and under the assumption of 
gaugino mass unification (black squares). 


of the gluino, and it is possible to build kinematic edges which can be used to constrain its 
mass. We can therefore assume that for the considered model it will be possible to exclude a 
significant stop coannihilation contribution to the relic density. The main loosely constrained 
parameters which influence the relic density calculation are therefore: tan/3, m{A), and m(f 2 ). 
We have therefore varied in turn m(f 2 ), tan/3, m{A), while the other soft parameters going in 
input to micrOMEGAs are recalculated in such a way that the masses of the sparticles which 
can be measured experimentally are kept fixed at their nominal value. 

For this exercise we have assumed the minimum constraint in the m{A) — tan /3 plane 
discussed above, assuming a lower bound on tan/3 of 3, and vetoing the regions where the 
heavy Higgs can be discovered by the LHC experiments at 5 ct through their SM decays. We 
have moreover imposed the loose bounds on m(f 2 ) discussed above. The dependencies on the 
different parameters are shown in Figures [TT)1 II 1 1 and II2L 

It is clear from these figures that if we can not extract from the data any direct information 
on the mass of the heavy Higgses, only an upper limit on the neutralino relic density can be 
given. If it is possible to exclude heavy Higgs masses lighter than 300 GeV through the peak 
search in the bb mass distribution, the spread in relic density measurement is of order 1%, and 
is strongly dependent on the experimental value of the lower limit on the mass of the heavy 
Higgs. If the heavy Higgses are directly observable through their decays to SUSY particles. 
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Figure 9: 


Value of the predicted relic density Vtydi^ as a function of the measured x? mass. 


the Higgs mass is fixed, and its experimental uncertainty does not contribute to the error on 
the relic density prediction, as for high Higgs masses the contribution of channels involving 
Higgses to the neutralino annihilation is negligible. 

In this case the dominant contribution to the uncertainty will come from the poorly con¬ 
strained value of tan (3. In the interval allowed by the non-observation of the SM decays of the 
Higgs, the relic density varies by ~11%, as shown in the right side of Figure ITTl It is interesting 
to disentangle the contributions of the different annihilation processes to the variation. In the 
right side of Figure we show the annihilation cross-sections (in units of their contribution 
to 1/H). For the different processes. The spread is dominated by the XiXi —> process. 

The reason is that for each value of tan (3 we recalculate the soft parameters in such a way that 
the sparticle masses, and the branching ratios which are measured experimentally are kept 
constant. Therefore the composition of the Xi and the value of Or vary as shown in Figures EJ 
producing the dependency observed in the full line in the left side of Figure 

An additional uncertainty will come from the value of m(f 2 ), which contributes a ~ 7% 
spread to the result, as shown in Figure II2L This is because the r 2 exchange contribution 
is opposite to the fi contribution. The cancellation appears in the s-wave part of the pair 
annihilation cross section, which is chirality suppressed. In the limit where the Higgsino 
component of the lightest neutralino can be ignored, the fj contribution to the amplitude can 
be expressed as 

M{swave) oc sincos[1/(1 + m^lrn^o) - 1/(1 -P m\/mf^o)]Zli (2) 
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Figure 10: Dependency of the relic density from the value of the pseudoscalar Higgs 

mass m{A). The dependency is shown respectively for the whole m{A) mass range (left), and 
for m{A) > 300 GeV (right). 


yielding the observed dependence of the annihilation cross-section on m{f 2 ). The mixing 
angle 6t- is kept constant by construction in Figure El because we use the ratio of decay 
widths F(x 2 ^ eex 5 )/r(x 2 ^ a-s the constraint on 0,-. 

5 Discussion 

In the above discussion we have considered a ‘bulk region’ SUSY model in which ff) annihila¬ 
tion in the early universe is dominated by diagrams involving light sleptons. This model lies 
within the mSUGRA sub-space of the full MSSM parameter space, however we have shown 
that it is possible to set useful constraints on the neutralino relic density without making 
special assumptions on the behaviour of the theory at high scale. 

It is interesting to consider whether this result applies more generally to SUSY models 
characterised by dominant light slepton annihilation. The starting point in the analysis 
is the possibility of isolating the decay chains ql —> X 2 Q (Jl —> X 2 d 

f^T^q T^T^q with sufficient statistics to be able to observe the kinematic edges providing 
the mass measurement. One therefore requires the mass hierarchy m{qL) > ni(x 2 ) > 
and equivalently for the fi. The former condition on m(x 2 ) a-iid is generally satisfied 

in the light slepton annihilation region, and if stau annihilation is to be relevant so must 
the latter. To avoid significant squark co-annihilation m{qi) must by definition be somewhat 
larger than the mass scale of the sleptons and lightest neutralino. It is therefore likely that at 
least the necessary decay chain(s) will occur, although the mass differences may be sufficiently 
small that one or more of the kinematic end-points are unobservable. 
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Figure 11: Left: Dependency of the relic density from the value of tan (5. Right: De¬ 

pendency of the annihilation cross-section for different processes, respectively XiXi 
(full line), XiXi (dashed line) and Xifi —> Z/At (dot-dashed line). The units are 

contributions to l/Qh'^. 


In order to obtain sufficient statistics of events containing the above decay chain, m{qif) 
must be < 1 TeV, and the y® ™ust have a large Wino component to couple to ql- Moreover, 
the mass of the left handed component of fi should not be too large, to avoid the decay 
X 2 ^ Fir saturating the X 2 branching ratio and thus killing the lepton signature. From the 
experimental point of view, we need to add the requirement that the masses of the involved 
sleptons are not too near to the mass of either neutralino. A further essential ingredient in the 
reconstruction of the neutralino mass matrix is the possibility of measuring the mass of the 
X 4 . This implies non-negligible gaugino components in the X 4 , the appropriate mass hierarchy 
and typically the knowledge of the £l mass, either through its direct production or through it 
appearance in a cascade decay IMI. Such requirements on the neutralino and squark sectors 
need not be satished in the light slepton Xi annihilation region; if this were indeed the case 
the analysis would be more challenging. 

Finally, we also require that it be possible to constrain alternative contributions to the 
annihilation cross-section through exclusion of the required values of tua and m{ti) by direct 
measurement or non-observation of a signal in an appropriate channel. If such constraints can 
not be obtained then the obtained value for Llh'^ is merely an upper limit. 

We may consider also whether the techniques described here may be applied to SUSY 
models in which alternative annihilation mechanisms dominate. In general in order to obtain 
an estimate of it is necessary not just to measure the dominant contributions to the 
X 4 annihilation cross-section but also to constrain all other possible contributions. In any 
such analysis therefore the masses and mixings of the neutralino, light slepton, stau, stop and 
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Figure 12: Dependency of the relie density from the value of the f '2 mass. 


Higgs sectors must all be measured or constrained. The goals of the analysis must therefore be 
similar to that described here and if accessible use will be made of the decay chains considered 
above. 

Although it is impossible to outline the required analyses for all possible models a few gen¬ 
eral observations may be made. Models in which stau co-annihilation dominates display similar 
phenomenology to the light slepton annihilation region, however the small mass difference be¬ 
tween the fi and x? makes observation and measurement of the qi X 2 Q r^r^q 

decay chain more difficult. In regions in which the x? possesses a significant Higgsino or Wino 
component (see e.g. EZIEH] ), leading to dominant (co-)annihilation to EW bosons the analysis 
would be qualitatively different to that described here and highly model dependent. If heavy 
Higgs annihilation is important techniques such as those described above for measuring its 
mass will be vital. When stop co-annihilation is enhanced the small mass difference between 
the strongly interacting ii and x? will likely make measurement of the masses very difficult 
indeed. 

6 Conclusions 

We have explored the potential of the LHC experiments for predicting the cosmological relic 
density of the LSP from detailed measurements of the SUSY spectrum and decay modes. No 
unification condition is imposed on the sparticle spectrum. We have focused on a model with 
essentially bino LSP and light sleptons, for which a detailed experimental study exists in the 
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literature. We have examined the relative roles of the different measurement uncertainties, 
and studied the uncertainties due to SUSY parameters which are relevant for the relic density 
calculation, and are poorly, if at all, constrained by the measurements at the LHC. For the 
experimentally accessible measurements, we have highlighted the role of the measurement of 
the rr edge. For the influence of badly measured parameters, the key issue is the observability 
of the heavy Higgses at the LHC. We have discussed in particular the prospect for observing 
the H/A in the cascade decays of the sparticles or in a SUSY decay mode. In case the position 
of the rr edge can be controlled at the level of ~1 GeV, and the LHC experiments can 
demonstrate that the Heavy Higgses H, A have a mass in excess of 300 GeV, the value of the 
neutralino relic density can be predicted as: 

= 0.108 ± 0.0l{stat + sys);[|;[j[] 2 (^(^))taoii(tan/ 3 );o;oo 5 (m(f 2 )) 

after three years of data taking at high luminosity, corresponding to an integrated luminosity 
of 300 fb“^. In case no experimental information on the heavy Higgs can be extracted from 
the LHC data, it will only be possible to put an upper limit of approximately 0.12 on the 
neutralino relic density. The discovery of the decay of the heavy Higgses to SUSY particles is 
however most probably statistics limited, and a luminosity upgrade of the LHC might allow 
the discovery of this decay mode, thus yielding also a constraint on tan (3. 
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